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ABSTRACT 



Context. Balmer lines serve as important indicators of stellar effective temperatures in late-type stellar spectra. One of their modelling 
uncertainties is the influence of convective flows on their shape. 

Aims. We aim to characterize the influence of convection on the wings of Balmer lines. 

Methods. We perform a differential comparison of synthetic Balmer line profiles obtained from 3D hydrodynamical model atmo- 
spheres and ID hydrostatic standard ones. The model parameters are appropriate for F,G,K dwarf and subgiant stars of metallicity 
ranging from solar to 10"^^ solar. 

Results. The shape of the Balmer lines predicted by 3D models can never be exactly reproduced by a ID model, irrespective of its 
effective temperature. We introduce the concept of a 3D temperature correction, as the effective temperature difference between a 3D 
model and a ID model which provides the closest match to the 3D profile. The temperature correction is different for the different 
members of the Balmer series and depends on the adopted mixing-length parameter auu: in the ID model. Among the investigated 
models, the 3D correction ranges from -300 K to +300 K. Horizontal temperature fluctuations tend to reduce the 3D correction. 
Conclusions. Accurate effective temperatures cannot be derived from the wings of Balmer lines, unless the effects of convection are 
properly accounted for The 3D models offer a physically well justified way of doing so. The use of ID models treating convection 
with the mixing-length theory do not appear to be suitable for this purpose. In particular, there are indications that it is not possible 
to determine a single value of (Imlt which will optimally reproduce the Balmer lines for any choice of atmospheric parameters. The 
investigation of a more extended grid and direct comparison with observed Balmer profiles will be carried out in the near future. 
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1. Introduction 

Balmer lines are prominent features in stellar spectra. There has 
been a long tradition of using the Balmer line wings for the deter- 
mination of eff ective tem peratures of F-K stars (a non-exhaustive 
hst includes: ICavreld e Strob ej 1 19601: Searle & Oke 1962; 
IStrohbachI 1 1 9701; iGehrenI 1 1 98 It ICavrel et al.i 1 1985i: iSoderblomI 
il986). However, the use of Balmer lines as temperature indica- 
tors requires fairly high quality spectra and a sophisticated theo- 
retical framework, both from the point of view of micro-physics 
and of the model atmo spheres employe d . It w as not until the 
comprehensive work of iFuhrmann et al] (11993) that the wings 
of Balmer lines became the temperature indicator of choice in 
many investigations. One of the advantages of Balmer line tem- 
peratures is that, unlike those based on colours or on the infrared 
flux method, they are reddening independent and may, in princi- 
ple, provide an accuracy of the order of 50 K. 

In the framework of ID homogeneous model atmospheres, 
in which convective ene rgy transport is describ e d by the mixing- 
length approximation (iBohm-Vitensd Il958l) . iFuhrmann et al.l 
(Il993h studied in detail the effects of convection on the wings 
of the Balmer lines. They convincingly demonstrated that, as a 
consequence of their different depth of formation, the response 
to a change in the adopted mixing-length parameter (Imlt of the 
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various members of the Balmer sequence is different. The first 
line of the series, H„, is fairly insensitive to the choice of ckmlt, 
at least for solar metallicity models, but the higher members are 
fairly sensitive. This provides, in principle, a way to select the 
most appropriate value of ccmlt- For a star for which the effec- 
tive temperature is known, like the Sun, one may select the cmlt 
which best reprod uces all member s of the Balmer series. 

At the time Fuhrmann et alj d 19931) performed their in- 
vestigation, hydrodynamical simulations in which convection 
is treated in a phys ical and non- parameterized way ha d 
just become available ( iNordlund & Stein 1991; Stefen|U991|). 
However, there were not enough simulations to cover satisfacto- 
rily the range of T^ff and log g which was pertinent to their inves- 
tigation; moreover those simulations were not very sophisticated 
in the treatment of line opacity. After 15 years the situation has 
greatly improved. Full three dimensional hydrodynamical simu- 
lations (hereafter 3D models, for short) and the associated line- 
formation codes have reached a level of sophistication in radia- 
tive transfer comparable to that of state-of-the art ID models and 
line-formation codes. More importantly, for the first time a fairly 
large grid of 3D models is available to allow a systematic inves- 
tigation of convection effects. 

To clarify how the results of 'Fuhrmann et al.J ( Il993h may be 
affected by inhomogeneities, we investigate Balmer hne forma- 
tion in hydrodynamical model atmospheres. Broadening due to 
convective velocities is only important in the inner line cores of 
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Table 1. Parameters of the 3D models. The grid resolution is set 
to 140 X 140 X 150 (A^x x A^y x N^), time indicates the total stellar 
time, and snaps the number of snapshots per model. 
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Balmer lines; in this sense, the outer parts of the Balmer lines 
are particularly well suited sensors of the thermal structure of an 
atmosphere. We use a sample of 3D model atmospheres cover- 
ing the range [M/H] = -3.0 to 0.0 in metalHcity, 5500-6550 K 
in Jeff ™d 3.5-4.5 in \ogg. Our investigation is purely theoreti- 
cal, and we express our results in terms of effective temperature 
differences between 3D and ID models. 

2. Model atmospheres & line formation 

The 3D model at mospheres we use have been computed with the 
CO^B OLD code dFrevtag et alj l2002l 120081: IWedemever eTal] 
12004 ). The code solves the time-dependent equations of com- 
pressible hydrodynamics coupled to radiative transfer in a con- 
stant gravity field in a Cartesian computational domain which is 
representative of a volume located at the stellar surface. Further 
d etails on computation al methods and validity tests can be found 
in'Ludw ig et al.l (l2009l) . In Table[T]we provide some of the basic 
properties of the 3D models employed in this paper 

In order to make our comparison strictly differential we em- 
ployed as a reference ID models computed with the LHD code, 
which employs the same microphysics and numerical scheme for 
radiative transfer as CO^BOLD. Further details on these mod- 
els can be found in Caffau & Ludwig (2007). These models are 
classical hydrostatic ID models; n o velocity fields are consid- 
ered, convection is treated using the l]VIihalasl(ll978h formulation 
of the mixing-length theory, and one has to assume a microtur- 
bulence parameter in associated spectral synthesis calculations. 
We refer to these models as IDlhd- 

In addition we employed the ID structure which is ob- 
tained by a temporal and spatial average of the 3D models 
over surfaces of constant Rosseland optical depth. We refer to 
these models as (3D). The (3D) model has, by construction, 
the mean temperature profile of the 3D model, but no horizon- 
tal temperature fluctuations are present. Therefore, a comparison 
3D - (3D) highlights the effects of such fluctuations, a compar- 
ison (3D) - IDlhd highlights the effect of differences in the 
temperature profiles, and the 3D - IDlhd comparison provides 
insight into the combined effects of different temperature pro- 
files and horizontal temperature fluctuations. 

The Balmer line profiles have been computed with the 
LinforSD cod^B In the version of the code use d here, the Balmer 
line profiles are computed with the theory of ICavrel & TravingI 
(11960) . In the future, more up-dated theories will be imple- 
mented in LinforSD. However, for the purpose of our strictly 
differential analysis, the theory employed should not be rele- 
vant. LinforSD is capable of computing line profiles both for 3D 
and ID models. The 3D synthesis is very CPU time demand- 
ing. For this reason, we computed for each line only a range of 

' fittp://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 
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Fig. 1. Hq, Balmer line profiles computed using a 3D model and 
ID models with varying mixing-length parameter a, for Teft = 
5920 K, log g = 4.50 and [Fe/H] = -3.00. 

5.5 nm around the line center at a resolution of 1.5 x 10"^ nm. 
Furthermore, we did not include any blending lines in the wings 
of the Balmer profiles. While this clearly limits the accuracy with 
which the higher members of the series can reproduce the obser- 
vations of more metal-rich stars, it is not relevant for the present 
differential theoretical comparison. 

3. 3D temperature correction 

To quantify the granulation effects, we use the concept of a 3D 
temperature correction. This is defined as the difference between 
the effective temperature of a 3D model and the temperature 
derived by fitting Balmer line profiles computed from a T^s 
sequence of IDlhd models with the same surface gravity and 
metallicity to the corresponding 3D profile. The 3D temperature 
correction is in general different for the different members of 
the Balmer series (here we investigate only the first three mem- 
bers) and depends on the (Tmlt adopted for the IDlhd models. 
Although we believe that the 3D temperature correction is a use- 
ful description, one should bear in mind that it is a simplification. 
In general, the whole profile computed from a 3D model has a 
different shape from that of a ID model; an example is given in 
FiglH No profile computed from a ID model, whatever the tem- 
perature, can exactly reproduce the profile computed from the 
3D model. The temperature correction singles out the ID model 
that provides the profile nearest to the profile computed from the 
3D model. Clearly, the concept of distance (near or far) has to be 
defined by a suitable metric. 

As a measure of the similarity, we used the root-mean-square 
deviation A^^^^^ between the normalized line profiles above a pre- 
scribed residual flux level, 

/-I 

where is the number data points making up the line profile, 
/^'-* the flux based on the 3D, and /"^ the flux based on the ID 
model. We added explicitly the Teff-dependence in the ID case 
as a reminder that we varied the T^s of the ID models in order 
to match the 3D result. The 3D correction is obtained as the dif- 
ference (in the sense 3D- ID) between the effective temperature 
of the 3D model to the best matching ID one. 
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Table 2. 3D corrections, AT, the quality of fit, QF, and the sensitivity of the fit, crj^ff for the different models in units of K. 
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To calculate A^^^^^ one has to define the portion of the profile 
over which this is computed. The line core must always be ex- 
cluded in these comparisons, since in real stars it is affected by 
the presence of a chromosphere. When fitting theoretical to ob- 
served spectra, the usual choice is to define a wavelength interval 
which defines the wings of t he line. Typical ch oices are 0.3 nm 
and 0.5 nm from line centre dCavrel et al.lll985l) . In the compar- 
ison of theoretical spectra, however, we may take advantage of 
the fact that the continuum is perfectly defined and consider the 
wings as the portion of the profile that is above a given thresh- 
old. In this way a similar fraction of the line wing is considered, 
irrespective of the temperature of the star. We believe this choice 
is the most appropriate for the comparison of theoretical spectra, 
and use a threshold of 0.8 in residual flux. 



In addition to the best matching T^s of a ID model, the fit- 
ting provided further information: i) the overall quality of the fit 
is related to the residual A^^^ at the best matching 7^^; ii) the 
sensitivity to which we can determine Teff is related to the rate 
of change of A^^, (more precisely its curvature) with respect to 
Jeff at the matching point. Since we are dealing here with the- 
oretical, essentially noise-less data, the sensitivity is not a real 
issue, and we can fix the best matching Teff to arbitrary preci- 
sion. However, in practice one is usually working with spectra 
of only finite signal-to-noise ratio (S/N) so that the sensitivity 
is important for the precision to which Teff can be obtained. To 
make a connection to the practical limitations, we present sensi- 
tivities and qualities of the fits in a form which makes it easy to 
relate them to the situation one encounters in finite S/N spectra. 



We assume a simple noise model where all pixels have the 
same S/N distributed according to Gaussian statistics. This is 
well justified since all pixels have a rather similar flux level due 
to our chosen high flux threshold in the fitting. The standard de- 
viation of the Gaussian distribution of the flux crp is given by 
CTp = (S /Ny\ Our minimization of A^^^ would then transform 
to a -minimization.;^^ is related to A^^^ i^sx^ - ^^rms where 

o-p 

is the number of points which sample the line profile. The un- 



certainty of the fi tted effective temperature crxeff is given by (e.g., 
iPress et al.|[T992l) 



0"Teff 



1 d^x 



2,2\ 



^' N 



2 \ 
Tins 



(2) 



Not surprisingly, the attainable precision scales inversely pro- 
portionally to the S/N, and to V/V. At this point, we perform a 
more subtle re-interpretation of A^: t2^^ - as an average - does 
not sensitively depend on A^ as long as the line profile is suffi- 
ciently densely sampled. Hence, Eq. (|2]) does not only provide 
the scaling of crxeff for the particular fit in question but for any 
fit based on a number of A^ points. We use this property to give 
CTxeff always for the same nominal S/N of 100 and number of 
points A^ = 700 making all presented fits inter-comparable. In 
the present context, cr-Yds is a measure of the sensitivity - or 
rather insensitivity - of the Balmer line profiles to changes in 
Teft . If one is dealing with spectra of finite S/N, crxeff provides 
an estimate of the attainable precision, and Eq. (|2]i can be used 
to translate the presented values to the actually present S/N and 
number A^ of statistically independent elements (pixels) in the 
spectrum under consideration. 

In a similar vein, we express the quality of the fit QF related 
to the residual A^^^^ at the fitting point in terms of a temperature 
difference 



QF 



, (residual) 



i2a2 \—i 



d'A, 



eff 



(3) 



This provides a handy measure of whether the derived 3D Teft 
corrections make sense; if the fit quality in terms of a tempera- 
ture difference QF is larger than the derived correction AT one 
should take this as a warning that the 3D versus ID differences 
of the line profiles are so large that trying to match them is not 
appropriate. 

4. Results 

For each 3D model, the 3D temperature correction (AT"), the 
quality of fit (QF) and the sensitivity (crxeff) for three different 
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Table 3. Spread of AT among Balmer series members. 
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grids of LHD models corresponding to three choices of omlt 
are listed in Table |2] The last nine columns of the table provide 
the same quantities, but for the (3D) model. These allow us to 
understand what part of the 3D temperature correction can be 
ascribed to differences in the mean temperature profiles of a 3D 
model and IDlhd model. Graphical representations of the re- 
sults are provided in the online appendix. 

The sample of investigated models is not large, however 
several patterns in the temperature corrections clearly emerge. 
At solar metallicity, aMLT= 0.5 provides nearly the same tem- 
perature correction for the first three Balmer lines. This re- 
sult matches what was derived from the comparison of ob- 
served sp ectra to theo retical spectra based on ID models by 
[puhrma nn et ai.l (119931) and Ivan't Veer-Menneret & Megessied 
tl996.) . We also find that the wings of Ha are relatively insensi- 
tive to ttMLT- In the past this has motivated the use of only Ha, 
but not of higher members of the series, for the determination of 
Teff ( Asplund et al. 2006; Bonifacio et al 2007). At low metal- 
licity, however, the insensitivity is lost, and the temperature cor- 
rection of Ho- is not independent of q-mlt- There are also sizeable 
differences in the temperature corrections of the various Balmer 
lines. 

In Table |3] we provide the temperature correction spread 6T, 
defined as AT^ax - AT^in among the first three lines of the 
Balmer series. For the models at reff=5500K and T^ff- 5920 K 
the smallest spread is achieved for q'mlt=2.0. For the T^^ff- 
5500 K model the spread in the temperature correction is not 
even a mono tonic function of ckmlt- However, the spread for 
Q'MLT=2.0 is not considerably smaller than that for Q'mlt= 0.5, 
especially when compared to the associated errors. This is not 
surprising; the mixing-length theory is a parametric phenomeno- 
logical description of convection and it is obvious that it is not 
possible for a single value of the free parameter to capture all the 
complexity of this physical phenomenon. It can be seen that the 
situation is more complicated for the cooler models, while for 
the three models hotter than 6000 K, 6T is a monotonic value of 
ajviLT and achieves the smallest value for aMLT= 0.5. 

The (3D) - IDlhd differences appear to be quite regular: 
the lower the temperature, the higher the temperature correc- 
tion, reaching the value of 500 K for the model at Teff =5500 K. 
This regularity is not immediately obvious when looking at the 
3D - IDlhd differences, since these depend also on the temper- 
ature fluctuation. The importance of temperature fluctuations is 
different for the different T^ff, but also for different metallicities 
and gravities. In all cases we see that the role of the temperature 
fluctuations is to reduce the temperature correction, with respect 
to what is derived with respect to the (3D) model. But this re- 
duction can be as large as 50% (for the Teff =5500 K model) or 
as small as 7% (for the reff=5920K model). 



5. Conclusions 

The main conclusion that can be drawn from this investigation 
is that stellar granulation has sizeable effects on the wings of 
the Balmer lines. When quantified in terms of the temperature 
correction, such differences span the range of +300 K for the in- 
vestigated models. This implies that if high accuracy effective 
temperatures are to be derived from the wings of Balmer lines, 
the effects of granulation must be taken into account. A temper- 
ature scale based on fitting the wings of Balmer lines with ID 
models will be different from that derived by using 3D models. 
The difference between the two scales is not a simple offset, but 
rather has a temperature dependent slope. 

The smallest spread in temperature correction often occurs 
when ffMLT - 0.5, but not always. The temperature correction 
for H„ is generally larger than that of the other members of the 
series. At lower temperatures the sensitivity of H„ to effective 
temperature drops. This is illustrated by our cooler model ij^s- 
5500 K); remarkably, the higher members of the series maintain 
a rather good sensitivity to T^g. The use of 3D models, removing 
the uncertainty on the choice of aMLT, suggests that the use of 
several Balmer lines should greatly increase the accuracy and 
robustness of the T^^ determination - provided one can handle 
the increasing line blending for the higher series members. 

While a bias of the true T^^ of a late-type star determined 
from ID models seems unavoidable, the question occurs as to 
whether a T^s corrected for the 3D- ID difference is really su- 
perior if one wants to use the underlying ID model for the in- 
terpretation of other features in the stellar spectrum - typically 
spectral lines for abundance determinations. Arguably, the cali- 
bration of the ID model inherent to the temperature fitting can 
make it advantageous to rather use the uncorrected ID temper- 
ature. However, this hinges on the relation between the spectral 
feature of interest and the one fitted for the Teff determination, 
and has to be decided upon on a case by case basis. In the on- 
line appendix we present a case where a Teff -corrected ID model 
performs better 

The general trends highlighted in our investigation need to 
be confirmed by the use of a larger grid of 3D models. In ad- 
dition the implementation of a more up-to-date line broadening 
theory will allow a direct comparison between 3D synthesis and 
observed spectra. In the near future we plan to extend our inves- 
tigation in these two directions. 
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Fig. A.l. The temperature differences (3D - ID) using ID models with a - 0.5, a = 1.0 and a - 2.0, for Ha, H/8, and Hy for all 
of the six models considered. The uncertainties related to each temperature difference measurement are displayed as boxes. The 
value of the uncertainty on the x-axis corresponds to the associated crjgff (sensitivity) values, while the uncertainty on the y-axis 
corresponds to the associated QF (quality of fit) values. See Table|2]and the text for details. 



Appendix A: Graphical representation of the 3D temperature correction, AT, the quality of the fit, QF, and 
the sensitivity of the fit, o-jeff 
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Fig. B.l. 3D-1D abundance correction versus effective temper- 
ature difference for the test case of a metal-poor F-type dwarf. 
For details see text. 

Appendix B: 3D-1D temperature correction and the 
lithium abundance in metal-poor F-type dwarfs 

Here we present an example where it is advantageous to use a ID 
model corrected for the 3D- ID temperature difference when de- 
riving the chemical abundance from spectral line analysis: the 
abundance of lithium in a metal-poor F-type dwarf obtained 
from the 670.7 nm resonance line. We performed a 3D-NLTE 
spectrum synthesis calculation for the line on the 3D model 
with reff=6280K, log 0=4.0, [M/H]=-2.0 (cf. Table[T]i. We con- 
sidered the resulting spectrum as representing an observation. 
Unlike a real observation, however, the underlying lithium abun- 
dance and stellar parameters are exactly known. According to 
Table [3] the temperature correction from Ha fitting amounts to 
74 K in this case, i.e. the ID model that fits the 3D Ho' pro- 
file best is 74 K cooler than the 3D model. We then calculated 
for a series of thirteen ID (LHD) models of different effective 
temperatures ID lithium line profiles in LTE and NLTE. In the 
ID spectrum synthesis, we assumed a microturbulence velocity 
of 1 km/s; however, the actual value is not important since the 
line was chosen to be very weak. We derived for each model the 
lithium abundance matching the line strength obtained in 3D. 
Figure [BTI depicts the resulting abundance differences between 
the underlying lithium abundance assumed in the 3D model and 
the derived ID abundance, versus the effective temperature dif- 
ferences between ID models and the 3D model. As evident from 
the plot, one reduces the abundance error resulting from the 
erroneous effective temperature of the ID model when apply- 
ing the 3D- ID temperature correction. This holds irrespective 
of whether the ID abundance analysis is performed in LTE or 
NLTE. Not surprisingly, the figure also shows that the correction 
of the effective temperature does not result in a perfect match of 
the lithium abundances in a ID analysis. 



